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. Strikes (RANS) equations are used to numerically simulate 

The three-dimensional, Reynokls-averag , ancle of attack. Two large-amplitude, high- 

nonacad) .ortic.l n»w about a ar , -u^aSld. The free-to-roti motion is eomputed 

the dvnamic forces and moments. 


Introduction 

T HE high angle-of-attack flight regime often includes 
complex phenomena such as nonsteady flow, crossflow 
separation, and vortex breakdown. Modern tacPcal fighters 
fly at high angles of attack in order to take advantage of the 
nonlinear lift generated from vortices that form on their le ^ 
ward sides. This results in a substantial improvement of a 
aircraft's maneuver and agility performance. However at suT 
ficiently high angles of attack, vortex asymmetries can form 
and induce dynamic motions such as wing rock a sustain 
limit-cycle roll and yaw oscillation. Wing rock and other self- 
induced aircraft motions can be difficult to control and may 
result in departure from controlled flight. 

A number of experiments 112 have investigated wing r 
and other dynamic motions for delta-wing geometries with a 
smgle-degree-of-freedom in roll. These simple geometries 
contain the relevant flow physics, e g., primal?, secondary, 
and tertiary vortices, and vortex breakdown without the ad- 
ditional complexities of a complete aircra ge 

perimental data usually include forces and moment^ as wel 

as on-surface and off-surface flow visualization. More re 
cently, some experiments 6 ' 2 have also measured time-de- 

^Cornp u tationa ^ fhi id 'd v n am ics (CFD) is just beginning to 
be used to investigate wing rock. CFD can provi e 
flowfietd detail than is otherwise possible by experiment a one, 
and therefore, complements experimental investigations. 
However CFD flow simulations can be computationally costly 
licause the flowfield is often nonsteady and highly nonlinear 
Some investigators have reduced the computational cost by 

assumption, similarity principles are used to reduce a three- 
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limensional physical flow into a tw0 ' di “rvortex^Sk- 
ion As a result, this approach cannot predict vortex breaK 
lown that often accompanies wing rock. Furthermore 
nethods based on the assumption of inviscid flow only r predtc, 
,ortex formation emanating from flow 
eading edges. The inviscid assumption completely ignores th 
/ortices that are known to form from separation at smooth 
surfaces e g , forebody, secondary, and tertiary vortices. Th 
strength'of the inviscid primary vortices is also Actionable 
and may affect the position of vortex breakdown. Ultima y, 
these computational simplifications will have on y res ric 
utility for complete aircraft applications. ^ • 

The approach adopted in this work is to utilize the Navier 
Stokes simulation (NSS) code to numerically compute h.gh- 
incidence vortical flow about a delta wing by solving the time- 

dependent three-dimensional, Reynolds-averaged. . avier 

Stokes (RANS) equations. Although this approach .s more 
computationally costly, it contains all of the relevant tlow 
P hS it is also anticipated that future computer improve- 
ments will continue to reduce the computer time needed to 
obtain time-dependent RANS solutions. 

The overall goal of this effort is to develop an expenmen- 
tally validated CFD tool for predicting and an C'? 8 * 
roll forced dynamic-roll motions, and free-to-roll motions for 
delta wings at high incidence using the RANS equations. 

vortical flow over a 65-deg sweep delta wing at °^ 7 
number and 15-deg angle of attack _The computed staPc and 
dynamic forces and moments, and time-dependent surface 
pressures were in very good agreement with large-amplitude, 

E!ra r S- m ouo„ extern,'” Thb ro,,;ng = c ffa 



dynamic case. 
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Fig. 2 Comparison of mean computed and experimental rolline-mo- 
ment Meffiamts for different static roll angles from Ref. 18. M. = 
U.Z7, a - 30 deg, Re = 3.67 x 10 6 . 

corresponds to an angular rotation about the wing's longi- 
tudinal axis. Note the very good agreement for static C, up 
through 42-deg of roll. The negative slope indicates this wing 
is stattcally stable in roll under the present flow conditions 

is a com 3 tnm P T' at * = 0 de §- Also shown in the figure 
is a comparison of computed and experimental C, for a large 

tinn P , H de , -“ 4 u deg) ’ h 'gh-rate (7-Hz), forced-roll mo- 

me eSc ts er ThVd Q VS * CUrve is due t0 roll 

te effects. The dynamic C, traces out the curve in a coun- 

erclockwise direction with increasing time, indicating that 

his wmg ,s positively damped in roll. The areas enclosed by 

the dynamic curves indicative of the aerodynamic damping 

coefficient matcheci within 3%. At this angle of attack no 

ei "“ “•”*««' or 

Chaderpanand SchifP- later computed the flowfield about 
the same 65-deg sweep delta wing at 30-deg angle of attack 
and sev , fixed roll angles. where vortex bre/kdoindoes 

^rh rnlT h ™ ean . stat,c rolling-moment coefficient variation 
with roll ang e is shown in Fig. 2. Note the nonlinear behavior 
of C, for small roll angles. This was attributed to vortex break 
l n h| T ree statica, 'y stable trim points and tSSSSt 
unstable tnm points are observed in both the computed and 
experimental results. The overall agreement is goodput some 
differences are apparent. It is anticipated that grid refinement 

Som 1 TMnnIl C ° mparis ° n ' Nevertheless, the present grid 
(about 700.000 grid points), was sufficient to predict the roll- 
ing moment nonlinearity. 0,1 

The purpose of this article is to build upon the previous 

tions'of the 65 j 0mpu, ' n S *be fl owfield for dynamic-roll mo- 
tions of the 65-deg sweep delta wing at 30-deg angle of attack 

where vortex breakdown does occur. Two forced-roll oscil- 
lations and a damped free-to-roll motion are presented The 

rions roThLfl SU h, S H re ° btained b - v coupling the RANS equa- 
tions to the flight dynamic equation of motion with a single- 
degree-of-freedom in roll. 8 

° f the wind ' tunnel experiment used to validate 
he computations is given later. This is followed bv a descrip- 
tion of the numerical approach, which includes the governing 
equations, turbulence model, numerical algorithm 8 comp U 8 
tational grids, and numerical boundary conditions. Compu- 

romarks * "* ^ presented and (ollow ed by concluding 

Summary of Wind-Tunnel Experiment 

Extensive wind-tunnel investigations for a 65-deo leadinp 

Hanff and his colleagues''-" in the Canadian Institute for 
Ae ™‘~ spaCe * e f arch 6 x 9 ft Low Speed Wind Tunnel and 
m the Wright Laboratory 7 x 10 ft Subsonic Aerodynamic 
Research Laboratory (SARL) tunnel. Very high-quality ex- 
perimental data was obtained for static (fixed) roll forced 
dynamic-roll motions, and free-to-roll motions in the form of 
mean and time-dependent surface pressures, forces, and mo- 



Fig. 3 Three-view drawing or the wind-tunnel model. 

"7** la ™ llght sheet technique was used to visualize 
nonsteady off-body vortices, whereas an oil flow technique 
was used to visualize surface-flow patterns for static roll cases 

These experimental data'" are used to validate the present 
flow computations. piesent 

F io A 3 th ™' V J 7 draW,ng ° f the de,ta -wing model is shown in 
*f_ : n T be d " lta u 7 n g ,s symmetric relative to its longitudinal 
axis and has double bevels that run parallel to and along the 
wing s leading and trailing edges. The delta wing is mounted 
a fuselage/sting that contains instrumentation and allows 
for a single degree of freedom in roll. The fuselage has a 
tangent-ogive cylinder shape and the sting has a reduced 

tathe flat^ ° f r 8 ! rail ' ng Cdge Beveis j° in the fuselage 
to the flat wing surface for additional stiffness. The model is 

constructed of a multilayer carbon-composite skin and foam 

moLnr r ng “■ 3 7 ff ' l, S htwei g ht mod el having a low rolling 
moment of inertia. Dynamic pressure transducers are located 
cross the span of the wing on the upper surface at 75% of 

found m Re°fs 6 9 rther det3i ' S 3b ° Ut ** eXperiment ca " be 

Numerical Approach 

Governing Equations 

The time-dependent RANS equations are transformed from 
Cartesian coordinates (x.y, z, and r) to body-fitted curvilinear 
coordinates (f ,, and r) to simplify the implementation of 
where 3 !? condltlo " sThe thm-layer approximation is used 

7hl hJ ‘ S aSSUr ? Cd that V,scous Ienms are m °st significant in 
the body-normal ^-coordinate direction. This approximation 
is consistent with the high Reynolds number flow condign" 
presented next The RANS equations can therefore be ex- 
pressed in the following strong conservation-law form: 


drQ + dfE + d v F + d,G = Re 


‘d c S 


( 1 ) 


p ,7 the '' ec,or of conserved dependent variables. E 
F and G are the inviscid flux vectors in the t v and > 
direcuons, respectively, and S is the thin-layer viscous flux 
vector. These equations have been nondimensionalized by the 
mg root chord C, the freestream density p„ and the free- 
stream speed of sound Any body movement or grid de- 
formation is accommodated through coordinate transforma- 
tion metrics that are included in the flux vectors. The perfect 
gas law, Sutherland s viscosity law, and a turbulence model 
completes the RANS system of equations. A description of 

deta k b f p nCe rn 0del IS g,ven in the next sect ' on - Further 
details ot Eq. (1) can be found in Ref. 19. 
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The NSS code currently has the capability to compute flows 
about a delta wing with a single degree of freedom in roll 
about its longitudinal axis. The roll angle for a forced periodic 
motion is given by 

<M0 = <t> 0 + tU* cos(wt) (2) 

where <!>„ is the mean or offset roll angle, < ^’! f t ! ie a „d Msthe 
of motion. o> is the circular frequency (<o - M). and 
time variable. The reduced frequency is defined by 


k = wbl2V„,_ 


(3) 


where b is the span of the delta wing and V* is the freestream 
Tied The roll angle for a free-to-roll motion, where the 
aerodynamics and body motion are coupled, is given by 


l<t> + C, sgn(d>) = KO 


(4) 

where / is the rolling moment of inertia, <*> = d<f>/dr and j 
= d 2 d>ldt 2 are the angular velocity and acceleration, respec- 
tively^ and /(/) is the instantaneous rolling moment obtain 
from the coupled time-dependent RANS equations. The sec- 
Si™ “dre left-hand side is a model fo, the mechantcal 

bearing friction, where 


sgn(d>) 



rja=SSSsSS2= 

Warming algorithm u - c ui oc k tridiagonal 

sipation and requires the so uuo^of a 5^^,^ ^ ^ in f plicit 

system of equations. The dl * g J cit seC ond/fourth-order 

dissipation * ha V s soluSon of ^alar pentadiagonal 
equations' 1 The diagonal algorithm 

gonthm and 1. it heret the diagonal algorithm 

However, Levy etal. naves attack js too 

advanced in time according to the expression 

= (2<t>" - <#>" ’) + (*' 2 ")['" " C ' s § n(Vr)1 (5a) 


Wonff 10 renorted that / = 0.27 lb-in. -s* - and Cj 4.0 lb in 
2, the bearing friction of the experimental roll mechanism 
I. be modeled as a small constant. Equations (2) and (4) 
define the roll angle for rigid body and gnd motions. 

Turbulence Model 

The Baldwin-Lomax algebraic modeW toget er wl 
Deeani-Schiff modification 21 " are used to model turbulent 
flow conditions. This is an efficient isotropic eddy-viscosity 
model that properly accounts for crossflow separation in the 
presenceof Strong vortices above the wing. The original un- 
modified Baldwin -Lomax model often chooses a 'eng* srale 
associated with a vortical flow structure that lies outside the 
viscous boundary layer. This results in an eddy viscosity that 
can be one or two orders of magnitude too large, and effec 
tivelv suppresses all but the most dominant vortex structure. 
The modified model of Degani and Schiff restricts the choice 
Jlength scalTto the boundary-layer region in a rational 
manner, and therefore, gives appropriate values for the eddy 

Vi Experimental surface-flow patterns 111 indicate that the 
Reynolds number based on the wing root chord (Re c - 3.67 
x V) is sufficiently high that turbulent transition »s confined 
to a very small region near the delta-wing apex. The com 
potations therefore assume fully turbulem flow^nning at 
the apex of the delta wing, i.e., transitional effects ig 

nored. 


where 


V<#>" = <t>" - <t>" ' 


(5b) 


Fauation (5b) is a first-order-accurate approximation to 
Equation (5a) is second-order-accurate in time when d> does 

not change sign and is 1 ’™“‘ as well 

The NSS code s spatial and time accuracy, 

as^zonal grid capability for treating complex geometries - 
has been previously demonstrated and reporte in e i 
ature. 

C TT^m“on., hyperbolic grid gener.ro,> « used 

axis extends upwind hofd 

mnr-chord length, in .he body-normal direcnon Thejnd 


N T^'i rTchote o. rwo implicit, approximately factored, 
cenmaldifference algorithms in rhe NSS code ro mtegmte Ekp 
m The choice of algorithm is case dependent. The first 
( 1 ) The ch f w j algorithm 2 '’ and the second is 

gonthm is the Beam warming _* algorithm due to 

l7, ag ° n tdXussee - Both Algorithms, as currently imple- 



pig. 4 Perspective view of the computational grid. 
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?Sh P °j ntS th , e ^ CirCUmferentlal dlrec,lon - and 49 points in 
the body-normal direction, totaling about 700,000 grid points 

There are more grid points on the leeward side of the wing 

to n r?so| C,r ?h n, | erent,a | d ' reCti0n) than on the windward side 
o resolve the leeward s.de vortices. This grid is identical to 

Figs Tand 2) W ' th NSS C ° de Refs ' 17 and 18 (see a,so 
The delta-wing grid was split into four zones in the stream- 
wise direction This was done in order to take full advantage 
f two available Cray C-90 supercomputers that had different 
amounts of main and secondary memories. 

Numerical Boundary Conditions 

"?' sl,p condition (zero velocity relative to a moving 
solid surface) , s imposed on the wing and fuselage/sting sur- 

tfon ThJ roi 3 ^ denS “ y 2nd pressure are f °und by extrapola- 
tion. The total energy per unit volume is then computed from 

die perfect gas law. Uniform flow is imposed at the far field 
7 a zer0 ' grad,ent condition is used at the outflow bound- 
ary. Flow variables are averaged across the wake cut that 
mends from the wing trailing edge downwind to the outflow 
boundary. Boundary conditions are imposed on the spherical 

fn tht y r aVera f " g n °r; ar,ableS (0ne grid P° int off *e axis) 
in the circumferential direction. ' 

Zonal boundary conditions are updated sequentially with 

he most recent data available. Since the grids have coincided 

surfeces with identical grid points, data transfer from one grid 

n ? “®. ther ,s accomplished by direct injection. A more com- 

£%££ST 27 0f ,he 2onal 


reduced frequency [Eq. (3)] of k = 0.20. There are 15 000 
time steps per cycle of motion. This corresponds to a non- 

ST. 0 00362 “ d * di ~*' 

Figure 5 shows the leeward-side vortices for the delta wing 
a L ' nStantaneoiiS ro)l an gles. The viewing angle if 

Thknn^ W ' n f K nd P er P end,cular to the zero roll planform 
portion of the roll motion corresponds to the second half 
of a periodic cycle Eq. (2)j where the right wing begins from 
r ~~ 28 - 2 de S and moves in a clockwise direction 

agam a7J ^ r °" and comes to rest 

, f ~ 28 2 de g- The maximum angular rate. <t> = 
1772 deg/s. occurs at zero roll angle. The vortices are nu- 
merically visualized by computing instantaneous streamlines 
emanating from the delta-wing apex. These streamlines cor- 
respond to a flow that is considered frozen at an instant of 

ticaUv deoict > ?h CeS V ' SUalized b ' V this approach should realis- 
wh tT h VO COreS upwind of vortex breakdown, 
where transverse movement of the cores is slow. Instanta- 

neous streamlines may be less realistic downwind of vortex 
breakdown, where the flow is highly nonsteady. The mstan- 
taneous location of the vortex breakdown positions and the 
corresponding direction of motion are indicated in the figure 
by arrows. Vortex breakdown is characterized by a rapid growth 
of the vortex core diameter and a rapid change in direct^ 
of the vortex core. Notice at <f> = - 28.2 deg (left wing down) 
vortex breakdown on the left wing occurs closer to fhe wing 
apex than does the right wing vortex breakdown location As 
the wing rotates clockwise, the left wing breakdown position 


Results 

Three dynamic-roll cases were simulated, two forced-roll 
oscillations and one damped frce-to-roll motion. usin° the 

SrjJSTS'J RANS eq a U f °" s The dire ««>" »' positive 

roll tor the 65-deg sweep delta wing is shown in Fig. 4. This 

corresponds to the wing moving clockwise (right wing down- 

ard) as viewed from behind the wing facing forward (pilot 

k In J C3 ^’ the freestream flow conditions were Mach 

number h' H V " 8 ' 6 ° fattack a = 30 d eg, a nd Reynolds 
number based on the wing root chord Re c = 3 67 x 10 6 

Experimental data provided by Hanff and Huang 1 " are used 
to validate the computed results. Experimental forces and 
momen s for the forced-roll cases are obtained by taking an 
ensemble average over several cycles of motion and then using 
a fast Fourier transform to represent the forces and moments 
in terms of the first 10 harmonics of the motion. Thus the 
time-history of the experimental forces and moments are a 
T " Cl ' 0n of the Fourier decomposition. The computed 
time-history results also use an ensemble average over four 
cycles of motion. No attempt is made to decompose the com- 
puted results into Fourier components. 

Each forced-roll case was started from a converged com- 
putation whose static roll angle corresponded to the maximum 
ro 1 angle of the oscillation, i.e., where the dynamic-roll os- 
c llation is momentarily at rest. This minimized the initial 
startup transients, which were found to completely decay dur- 
ing the first cycle of motion. Thus, the computations' were 
found to be periodic during the second and following cycles. 

Forced-Roll Oscillations 

= 0 Deg 

The aerodynamic response for a 65-deg sweep delta wing 
o nd t rg0,nga 'arge-amphtude, high-rate roll morion about it! 
longitudinal axis -s obtained by numerically integrating the 
RANS equations with the NSS code subject to a forced rigid 

oSt g a7l7r77 nedby E< J; (2) - The fir st case has a zero 
onset angle (<t>„ - 0 deg), a roll amplitude $ = 28 2 dep 

and a frequency/ = 10 Hz. This frequency coTresponds ft 
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Vortex Breakdown 




Fig. 6 Helicity -density visualization of crossflow vortices at ^ 

16 9 deg. M. = 0.27, « = 30 deg, *„ = 0 deg, “eg, 

Ic = 0^20, and Re = 3.67 x 10*. xlc = a) 0.479 and b) 0.506. 


first moves towards the wing apex, reverses direction near the 
aD ex and then proceeds downwind towards the wing trailing 
edge’ Conversely, the right wing breakdown position first 
moves downwind toward the wing trailing edge, reverses i 
rection near the trailing edge, and then moves upwind towards 
teling apex. The mo, ion of vertex breakdown ,s therefore 
seen to traverse a large portion of the delta-wing length. There 
is a clear lag in the vortex breakdown motion relative to the 
forced-roll motion of the delta-wing body. As expected, the 
vortex breakdown positions on each side of the wing move 
in opposite directions relative to each other. Figures 5d an 
5e show a large change in the instantaneous vortex breakdown 
positions, when the wing rolls from <f> - 12.6 to 28.2 deg 
This indicates the vortex breakdown positions move at a very 
high rate during this portion of the roll oscillation. The in- 
stantaneous streamlines in Fig. 5 imply spiral vortex brea - 
down. A separate nonsteady numerical flow visualization of 
the vortex core confirmed the presence of spiral vortex break- 
down This flow visualization was carried out in a postpro- 
cessing manner by releasing particles in the vortex core, near 
the delta-wing apex, and tracking their motion in ime by 
integrating the nonsteady Navier-Stokes velocity field. 

Figure 6 shows the instantaneous structure of the vortices 
in crossflow planes located just upwind and downwind of the 
right wing vortex breakdown, when the instantaneous roll 
angle is 16.9 deg and the wing is rotating clockwise. The 
streamwise location is approximately at 50% of the root chord. 
Helicity density-'' is used to visualize the vortices and is de- 
fined by 


H, = V ft 


(6) 


This is evident by the small pocket of reveled ^“ hed 
contours) in the right-wing vortex core. The sparseness oi 
helicity-density contours at the right-wing vortex core also 

indicate its weakened state. t 

The computed and experimental 10 dynamic rolling-momen 
coefficients C, are shown in Fig. 7. Computed static , va ^ues 
from Ref. 18 are also shown in Fig. 7, with additional st 
cases at d> = ± 65.6 deg. Both the computed and experimental 
dvnam c C are in good agreement with each other. Reducing 
IS SjSy . tow of 2 had no significant effect on the 
dynamic rolling-moment time history. Notice that * 

namic C, is not simply a perturbation superimposed on the 
stafic C, curve, as was the case in Fig. 1 where the flow had 
no vortex breakdown. The nonlinear time lags of the vortex 
breakdown position with respect to the delta-wing motion (see 
Fig 5) are responsible for this marked differen 
enclosed by a dynamic C, curve indicates the amount of non- 
dimensional work done by the fluid on .the • ™ng.The com- 
nutational and experimental nondimensional work are 3. 
and - 2 4 deg, respectively. Moreover, both dynamic C, curves 
trace out counterclockwise paths with increasing time as in- 
dicated by the arrows in Fig. 7. The negative work and coun- 
terclockwise paths indicate the wing is positively damped under 
the" present flow conditions. Since the work ***** - h 
the numerical flow is greater m magnitude than e s 
associated with the experimental flow, the computed flow 
more highly damped than the experiment. 

Figure 8 compares the computed and experimental dynami 
normal-force coefficients, C N . They are in fair agreement^ 
Both curves have the same shape, but the computed norm 
loads are less than those found in the e^ment. Thi 
consistent with the static results shown in Ref. 18 where 
computed suction peaks downwind of vortex breakdown were 

less Than those found in experiment. Utilization of a f,ner S" 
in the crossflow direction should improve the comparis . 
The computed and experimental 



Fig. 7 Dynamic 
a = 30 deg, <*> 0 = 
x 10 6 . 


static rolling-moment coefficients. M. - 0.27, 
ieB <t> = 28.2 deg, * = 0.20, and Re = 3.67 


where V is the fluid velocity and ft is the fluid vorticity. A 
strong vortex will tend to have high-valued contours clustered 
about the vortex core, whereas a weaker vortex will tend to 
have low-valued contours sparsely spaced about the vortex 
core. The sign of helicity density indicates th ® dlrect ,f" • 
vortex rotation, when there is no vortex breakdown. Th 
allows the easy identification of primary, secondary, and ter- 
tiary structures. However, the helicity density of a vortex core 
can change sign when vortex breakdown occurs^ This is due 
to the presence of local reversed flow. In Fig. 6, solid line 
indicate positive helicity density, while dashed lines indicate 
negative helicity density. Primary, secondary, and tertiary 
"offices are present in the flow. In Fig. 6a there ,s vortex 
breakdown above the left wing while the right-wing vortex is 
still intact. In Fig. 6b, which is located slightly downwind of 
Fig. 6a. vortex breakdown occurs on both sides of the wing. 



t Dynamic normal-force coefficients. M. 30 eg, 

0 deg, = 28-2 deg, * = «•“- and Re = 367 X 
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Fig. 9 Dynamic streamwise c.p. locations. = 0.27, a - 30 deg, Fig. 11 Dynamic normal-force coefficients. = 0.27, a = 30 deg, 
= 0 deg, = 28.2 deg, k = 0.20, and Re = 3.67 x W. <J> 0 = 28.0 deg, = 31.9 deg, A: = 0.20, and Re = 3.67 x 10 6 . 



Fig. 10 Dynamic and static rolling-moment coefficients. = 0.27, 
a = 30 deg, 4> 0 = 28.0 deg, = 31.9 deg, k = 0.20, and Re — 
3.67 x 10 6 . 



Fig. 12 Dynamic streamwise c.p. locations. = 0.27, cr = 30 deg, 
<D 0 = 28.0 deg, = 31.9 deg, k = 0.20, and Re = 3.67 x 10*. 


the computed X c p located slightly upwind of the experimen- 
tally measured values. 

<!>„ = 28 Deg 

A second forced dynamic-roll motion was computed with 
an offset angle <1> 0 = 28.0 deg, a roll amplitude 3> max = 31.9 
deg, and / = 10 Hz (k = 0.20). This corresponds to a roll- 
angle range of -3.9 ^ <f> ^ 59.9 deg. The nondimensional 
time step (A r = 0.00362) was identical to the previous case. 
The computed and experimental 10 dynamic C, are shown in 
Fig. 10, together with computed static C, values from Ref. 
18. Arrows in the figure indicate how the experimental C { 
vary with increasing time. Note that there is a large counter- 
clockwise loop and a small clockwise loop. The counterclock- 
wise loop corresponds to a damped portion of the wing mo- 
tion, while the small clockwise loop (55.2 < </> < 59.9 deg) 
corresponds to an undamped portion of the wing motion. 
Good agreement is seen between the computed and experi- 
mental C,. The dynamic C } curves differ even more dramat- 
ically from the static curve than in the previous case. This 
behavior is again attributed to the large time lags in vortex 
breakdown motion. The nondimensional work done by the 
fluid on the wing is again evaluated by computing the net 
area enclosed by the dynamic curve, with the area inside the 
clockwise loop considered as positive. The computed and ex- 
perimental net nondimensional work are -1.8 and -1.4 deg, 
respectively. These values, together with the dominant coun- 
terclockwise loop, indicate that this wing is also positively 
damped, as in the previous dynamic case. Moreover, the com- 
puted flow is again found to be more highly damped than the 
experiment. 

The computed and experimental dynamic C N are shown in 
Fig. 11 and are in good agreement with each other. Both 
curves form a figure-eight pattern. Note that the smaller C N 
values in this dynamic case result in better agreement between 
computation and experiment than in the previous dynamic 
case, which had larger C N values. This again suggests that 
further improvement could be obtained by refining the grid. 


Figure 12 shows the computed and experimental dynamic 
X c p , which are in good agreement with each other. 

Damped Free-to-Roll Motions 

A final dynamic case is presented in which the delta wing 
is initially fixed at <f> — 40.5 deg and then released. The wing 
is allowed to roll freely under the influence of the instanta- 
neous aerodynamic rolling moment. The time-dependent RANS 
equations are therefore coupled with the flight dynamic equa- 
tion of motion that has a single degree of freedom in roll, 
[see Eqs. (1) and (4)]. The details of this coupling were pre- 
viously described in the Numerical Approach section. The 
nondimensional time step was At = 0.005. 

Figure 13 shows the time-history of the computed and 
experimental 6 roll angles. Also shown in the figure is a roll 
simulation using a linearized model that was reported by Hanff. 6 
The computed and experimental roll angles are in close agree- 
ment with each other for the first half-cycle of roll motion. 
Thereafter, the computed roll motion damps out more rapidly 
than the experiment. This result was anticipated, since both 
of the computational forced-roll motions described earlier 
exhibited greater damping than the experimental motions. 
The locations of the minimum and maximum peaks in the 
computed roll-angle oscillations are nearly identical to those 
found in the experiment. This indicates the computed fre- 
quency of motion is close to the experimental frequency. On 
the other hand, the linearized simulation of Hanff 6 has a lower 
frequency and takes much longer to damp out. The final 
computed and experimental roll angles, where the delta wing 
comes to rest, were 17.5 and 20.3 deg, respectively. The com- 
putational static C ( curve in Fig. 2 has a zero rolling moment 
at <p ~ 16 deg. The difference between this static roll angle 
trim point and the final static roll angle in the damped free- 
to-roll computation is due to the friction term in Eq. (4). 

The phase-plane curve <£(/) vs <f > , is shown in Fig. 14. This 
trajectory begins at <p = 40.5 deg and <f> = 0. The curve 
follows a clockwise path until it reaches the final static con- 
dition at <f> = 17.5 deg. The final roll angle for the free-to- 
roll experiment is also shown in the figure. 
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Fig. 13 Time history of roll angle for free-to-roll motion. M„ = 0.27, 
a = 30 deg, <l> 0 = 40.5 deg, and Re = 3.67 x 10*. 



Fig. 14 Phase-plane path for a damped free-to-roll motion. Af. = 
0.27, a = 30 deg, <i> 0 = 40.5 deg, and Re = 3.67 x 10*. 

Conclusions 

The dynamic flowfield about a 65-deg sweep delta wing at 
30-deg angle of attack and Re = 3.67 x 10 6 was numerically 
simulated using the time-dependent, three-dimensional, RANS 
equations. Two large-amplitude, high-rate, forced-roll mo- 
tions, and a damped free-to-roll motion were investigated. 
The latter computation coupled the RANS equations to the 
flight dynamic equation of motion with a single degree of 
freedom in roll. Vortex breakdown was present in each case. 

The forced-roll computations were validated by comparing 
the computed and experimental normal-force coefficients, 
rolling-moment coefficients, and streamwise c.p. locations. 
The overall agreement was good, with some differences at- 
tributed to grid coarseness. There were large time lags in the 
vortex breakdown motion relative to the forced-roll motion 
of the delta-wing body. This resulted in a complex hysteresis 
of the dynamic rolling moment coefficient that differed sub- 
stantially from previously published static results. The wing 
motion was positively damped in both cases. 

A free-to-roll motion was also computed, where the delta 
wing was initially released from rest at 40.5 deg of roll. The 
resulting computational free-to-roll motion damped out more 
rapidly than the experiment, but the computed and experi- 
mental frequencies were in good agreement. There were three 
possible statically stable trim points. The computed and ex- 
perimental roll motions went to the same trim point. 
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